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We explore novel topological phases realized in a superlattice system based on the Wilson-Dirac model. Our main fo- 
cus is on a two-dimensional analogue of the weak topological insulator phases. We find such phases as are characterized 
by gapless edge states that are protected by symmetry but in a way sensitive to the relative orientation of the edge to the 
superlattice structure. We show that manifest and hidden reflection symmetries protect such weak topological phases, 
and propose bulk Z 2 indices responsible for the topological protection of the edge states. 
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Topological classification is a new trend in the field of con- 
densed matter. Though realized recently to the wide com- 
munity, its position is more than influential for determining 
future directions of the field. Beyond the Ginzburg-Landau 
paradigm, topological classification can be applicable for 
quantum liquids without fundamental symmetry breaking. 12 
Still truly generic systems are boring and the symmetry again 
restricts the systems. Then we have various physically inter- 
esting phases protected by symmetries. The periodic table of 
topological phases as an extension of the classical symmetry 
classes is situated at the heart of the idea. 3-8 Recently, exten- 
sion of the standard classification scheme by including the di- 
versity of topological phases protected by other types of sym- 
metry has been investigated. 9 ~ 16 

Topological phases are mostly gapped then the bulk is char- 
acterized by the absence of low energy excitations. On the 
other hand, with boundaries or impurities, there exist peculiar 
localized modes as edge states. The emergence of the edge 
states is not accidental and is a fundamental property of the 
topological phases as the bulk-edge correspondence. 7, 8 ' 17 The 
edge states further reflect symmetry of the topological phases 
and describe their variety beyond the bulk characterization. 

In this letter we attempt to further extend the idea of 
the topological classification to a superlattice version of the 
Wilson-Dirac type lattice model that exhibits hidden reflec- 
tion symmetry. Motivation of this work is by no means purely 
academic. Now that basic understanding of simple topolog- 
ical insulator crystals has been established, a possible direc- 
tion of not only the theoretical but the experimental research 
is to seek for the diversity of topological quantum phenom- 
ena. Recently, multilayer heterostructures consisting of alter- 
nating layers of topological and ordinary insulators have been 
experimentally realized, exhibiting an interesting correlation 
of bulk and surface properties. 18, 19 

This letter highlights a two-dimensional (2D) analogue of 
such a superlattice system; a variant of the Wilson-Dirac type 
tight-binding model with a stripe structure (see Fig. 1). With- 
out such a spatial non-uniformity the simple 2D Wilson-Dirac 
model is known as a typical Z-type model specified by the 
Chern number. 20 Extension of this model to the quantum spin 
Hall effect (QSHE) 21 has been carried out, 22 in which the Z 2 
invariant 23,24 distinguishes the QSHE phase and trivial phase. 
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Fig. 1. The superlattice structure: (a) schematic illustration of the lattice, 
(b-c) different types of ribbon geometry. The two edges of the ribbon are 
along the y-axis in (b), while along the x-axis in (c). 



Here, we demonstrate that the superlattice version exhibits a 
richer phase diagram (see Fig. 2) that cannot be classified by a 
single topological invariant. Since the mass parameter m con- 
trols the Chern number c in the uniform system, the superlat- 
tice of the mass (ra + ,ra_) can be regarded as that of distinct 
Chern insulators (c+,c-). On the other hand, this model has 
the total Chern number C in its own right. We find that in the 
C - sector there appear a 2D analogue of weak topological 
insulating phases 25-27 characterized by anisotropic topologi- 
cal properties. Although it is gapless as a two dimensional 
system, we mention here similarity of graphene where there 
exist direction dependent boundary states with time reversal 
symmetry (C = 0). 28 - 29 

Figure 1 shows schematic configuration of the superlattice 
model considered from now on. The model is defined on a 
square lattice, and on each site r - (x, y) of the lattice the 
electron is allowed to occupy two orbital (pseudo-spin) states 
with which a set of Pauli matrices 07, (ju - x, y, z) is associ- 
ated. The Hamiltonian is defined by 



H = Z Z ( |r)r " <r + fl + |r + ^ (r| ) + i r > y w< r i 



(1) 



where jl stands for the unit vector to the /i-direction, and the 
hopping and the on-site potential terms are specified as 

if b 
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Fig. 2. Phase diagram of the superlattice model for b = 2 and t = 1 : Chem 
number C evaluated at each point on the (m+,m_)-plane is indicated by dif- 
ferent colors: blue: C = 1, red: C = -1 and white: C = 0. Thick black dots 
correspond to the reference points in Table II. 



Fig. 3. Energy spectrum in the ribbon geometry: b = 2, t = 1 with bound- 
aries along the y- (left) and x- (right) directions. The upper [lower] panels: (a), 
(b) [(c), (d)] correspond to the case of mass parameters m + = 2.5, m_ = 5.5 
[m+ = 3.5, m_ = 7] belonging to the phase A [B]. In (a), the suffix L and 
R indicate that the corresponding edge state is localized at the left or right 
boundary in Fig. 1 (b). 



V(r) = [m + (-l) x 6m - 2b] cr : = (m ± - 2b) cr z . (2) 

Notice that the lattice constant has been chosen to be unity, 
and hence, x and y take only integral values. In this letter we 
focus on the simplest superlattice structure as depicted in Fig. 
1 in which V(r) takes two alternating values on each column 
of the vertical stripe. Taking into account a doubling of the 
unit cell due to the stripe texture as a sub-lattice degree of 
freedom, one can block-diagonalize eq. (1) in the reciprocal 
space as H — Y>k \k)'H{k){k\ whose matrix element 'H(k) is 
given by 



<H{k) 



M + + t sin kyCTy 
r| + e m *T x 



r, + e~ m *rl 

M_ + 1 sin kyCTy 



(3) 



where different rows and columns specify the sub-lattice, i.e., 
whether the electron is on a blue column (x: even) or on a red 
one (x: odd) in Fig. 1 (a), and M ± - \m ± + b(cosk y - 2) cr : . 

The uniform line, ra+ = m_ = m (6m - 0), corresponds to 
the standard Wilson-Dirac model: 20 ' 22 'Hik) - tzZ/i sm ^< cr /j + 
m(k)cr z , where m(k) — [m + £>2/.<( cos ^u ~~ l)]°"c- The half- 
filled ground state of the model is classified by the Chern 
number 30 ' 31 C that takes a nontrivial value C — 1 when 
< m/b < 2, while C - -1 when 2 < m/b < 4 [other- 
wise C takes a trivial value (C = 0)]; change of the topo- 
logical number corresponds to closing of the gap [zeros of 
m{k)] at the Dirac point k = (0, 0) and at its doublers' points 
(7T, 0), (n, n), (0, n). Away from the uniform line it is still pos- 
sible to compute the Chern number C by applying the pre- 
scription given in Ref. 32 to the present superlattice model. A 
phase diagram thus obtained is shown in Fig. 2. 

It has the following specific features: Topologically non- 
trivial phases with nonzero Chern numbers C = ±1 extend 
from the uniform line to a region of m+ + m_. The regions of 
C — ±1 overlap; at least they look like doing so in the phase 
diagram, to form a finite domain of C — phase represented 
by A in Fig. 2. There appear other C — phases in different 
parts of the phase diagram separated by topologically non- 
trivial phases. Are these C — phases simply topologically 
trivial? Our answer is "No", in the phases represented by "A" 
and "B". These phases exhibit gapless edge states in the rib- 
bon geometry [panels (a) and (d) of Fig. 3], Interestingly, the 
way these edge states appear depends on the way the system's 



boundaries are introduced (compare the top and bottom rows 
of Fig. 3). The structure of the phase diagram, especially the 
shape of the regions A and B is strongly dependent on the 
hopping amplitude t in contrast to the phase boundaries on 
the uniform line: m+ - m_. 33 It should be emphasized that 
though the concrete arrangement of distinct regions shown in 
Fig. 2 is not generic and varies continuously as a function of 
t, the behavior of the edge states in the two types of C — 
phases A and B are generic, implying that they are protected 
by some symmetry. 

The concept of the bulk-edge correspondence is now es- 
tablished, 17,34_36 dictating that the topological property of the 
bulk is fully reflected in the spectrum of the edge states. In the 
remainder of this letter we show that the gapless edge modes 
found in the C — phases A and B are indeed topologically 
protected, and interpret the corresponding C — phases as 
a 2D analogue of the weak topological insulator in 3D. 25-27 
Figure 3 highlights the "weak" nature of the phases A and B. 

In the phase A, edge states appear at the boundaries parallel 
to the y-axis [Fig. 3 (a)], whereas no edge states appear along 
those parallel to the x-axis [Fig. 3 (c)]. As indicated in Fig. 
3 (a) a pair of counter-propagating modes; one localized at 
the left boundary (L) and the other at the right boundary (R), 
cross at the zero energy, and also at a specific value of the 
momentum (symmetric point), either at k y - or k y = n. As 
we will show soon, this phase is reminiscent of the QSHE, 
since C - c + +c- — for c + — 1 and c_ = — 1, and one pair of 
edge states are due to c + = 1 , and the others due to c_ = - 1 . In 
passing, we mention that in the present model, the existence of 
particle-hole and inversion (or reflection) symmetries ensures 
the spectrum at each k to be symmetric with respect to zero. 
This is also the case with the edge spectrum. 

The phase B is, on the other hand, understood by consider- 
ing the limit m_ — > oo. We will claim that the mid-gap states 
in Fig. 3 (d) can be deformed into topologically protected flat 
bands in this limit. To be concrete, when m_ — » oo, electron 
occupation on m_ sites is suppressed, and the model reduces 
to just a set of isolated one-dimensional ladders described by 



J. Phys. Soc. Jpn. 



LETTERS 



Table I. Non-zero coefficients of the Hamiltonian (6). 
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P x <H(k X ,k y )Pl l =<H(k X ,-ky), 



(9) 



a reduced Hamiltonian, 

'Hiky) - t sin k y (T v + [m+ + b(cos k y - 2)1 cr z , (4) 

Owing to chiral symmetry of this Hamiltonian, the Berry 
phase integrated over k y is quantized to or 7r. 29 To see this, 
set Y - t sin k y and Z = m+ + b(cosk y - 2). Then, if the ori- 
gin (0,0) is located inside the ellipse (Y, Z) forms in the Y-Z 
plane, i.e., 1 < m + /b < 3, the Berry phase is n (nontrivial) in 
which the zero energy flat bands are expected. 29 Thus the iso- 
lated mid-gap states in Fig. 3 (d) can be deformed without the 
bulk gap closing to these topologically protected flat bands. 
Indeed, in the phase diagram in Fig. 2 the blue and red re- 
gions become narrower as m_ — > oo, converging respectively 
to a linear region on m + = 2 and on m + - 6. 

What is a bulk topological invariant characterizing these 
weak topological phases embedded in C — 0? In contrast to 
the case of the so-called Z2 topological insulator the present 
system lacks time-reversal symmetry. Yet, as we demonstrate 
below, the proposed weak topological phases are protected by 
another type of a Z2 invariant associated with manifest as well 
as hidden reflection symmetry. To see this, it is convenient to 
introduce a unitary-transformed Hamiltonian: 



<H(k) = U(k x )<H(kW(k x ), 

U(k x )= l 2 (g)diag(l,e- iA >), 



(5) 



where II2 operates on the Pauli matrices. The transformed 
Hamiltonian 'H(k) is represented simply as 



'H(k) = 0;^, +^72^, 



(6) 



where 7-matrices are defined by j\ - cr x ® cr x , 72 = cr y <8> I2, 
y 3 = <t z ®1 2 , 74 = o- A (g>cr v , y 5 = cr x ®cr z , andy A , v = i[y M ,y v ]/2. 
Coefficients a^ and yS A , are listed in Table I. 

We begin by demonstrating the following two properties: 
(i) < H{k) possesses not only particle-hole symmetry but also 
reflection (inversion) symmetry, and (ii) 'H{k) is not periodic 
with respect to k x , fi(k x + n, k y ) - U(n)1-((k x , kyjU^in). Here, 
the form of U(n) — 1 2 ® 0% implies that a twisted boundary 
condition is imposed on the Hamiltonian < H{k). Nevertheless, 
the half-filled ground states of < H{k) and 1-l{k) give the same 
Chern number on the same Brillouin zone [0, n] <g> [0, 2n~\. 

(i) Let us first note that the Hamiltonian (6) has particle- 
hole symmetry, 



E<H(k)Z- 



-'H(-k), 



(7) 



where S = -i7273^ = cr\ <8> ioK and K is the complex conju- 
gation operator. The presence of this symmetry is rather natu- 
ral if one recalls that eq. (1) is a straight-forward extension of 
the Wilson-Dirac Hamiltonian. The model has other symme- 
tries prescribed as 



P y ?i(k x ,k y )P y : 1 ='H(-k x ,k y ), 



(8) 



where P y - y^K and P x - K. These may be regarded as 
(anti-unitary) reflection symmetry with respect to y- and x- 
direction, respectively, and therefore, the model has inversion 
symmetry P<H(k)P- 1 = f((-k), where P = P x P y = 73. Sym- 
metries (7) and (8) are manifest also in 'Hik), but the symme- 
try (9) is hidden in Hik). 

(ii) We introduce the Chern number C in the transformed 
basis. Let \fi{k) be the negative energy multiplet of the Hamil- 
tonian n-i{k) with a phase convention, 



lff(-k X ,ky) = P y lk(k X ,ky), 
lf{k X ,-ky)=P x lf{k X ,k y ). 



(10) 



Note that the periodicity of the wave functions is such that 
t//(k x + n,k y ) — U(n)tf/(k x ,k y ), where U{n) - -\y\j\. Let 
A, = tfr^djiff and F12 = etfidiAj be the Berry connection and 
curvature, respectively, where <9, = dk r Because of eqs. (10), 
these obey 

Aj(-k x ,ky) = (-y-'Ajik^k,), 

Aj(k x ,-k y ) = (-yAj(k x ,k y ), 

F~n(-k x , ky) = F 12 (k x , -k y ) = F 12 (k x , k y ). (11) 

The Chern number C of the half-filled states is given by the 
integration of F^ik) over the Brillouin zone. Let us verify 
C — C, where C is defined in terms of the wave function if/(k) 
of the original Hamiltonian 1-i{k) in (3). The two wave func- 
tions can be related as 4>{k) - U{k x )i[>{k). This implies A,- = 
Ai + ifStfdiUifr, and hence, F u = F n + eijditrP.U^djU), 
where A, and F\ 2 are the Berry connection and curvature de- 
fined through ifrik), and P-(k) is the projection operator to the 
occupied states, P-(k) - t(r(k)i// T (k). Then, since U~diU = 
-U2 ® diag(0, 1) is a constant matrix and P-(k) is gauge- 
invariant as well as periodic on the Brillouin zone, the above 
difference term between F and F vanishes if it is integrated 
over the Brillouin zone owing to the Stokes theorem on the 
torus. Thus, we reach C - C. 

Based on the properties (i) and (ii), let us consider a topo- 
logical invariant that characterizes the weak topological phase 
studied so far. Consider the Berry connection A,. Because of 
(10), the obstruction of the gauge-fixing occurs mainly along 
the four symmetry lines k x - 0, n/2 and k y - 0, n. Namely, if 
the Berry connection is plotted on the Brillouin zone, vortices 
can appear on these lines. Moreover, they are always paired 
because of (11). For example, on the Brillouin zone defined 
by [0,7!"] ® [0, 2tt], a pair on k x = line sit on the points 
(0, £*) and (0, 2n - k*) with the same vorticity. Of course, 
in some other points away from these symmetry lines, there 
can also occur obstructions. In this case, vortices appear "in 
quartets" which are symmetric with respect to four symme- 
try lines. These have also the same vorticity. Thus we know 
that even number of vortices always appear as long as the 
phase convention (10) is adopted. However, there are excep- 
tions: The four crossing points of the four symmetry lines, i.e., 
(0, 0), (n/2, 0), {71 12, 7r), and (0, n), which will be referred to as 
Xi, ■ ■ ■ , X4 in this order. On these points, single vortices can 
appear. These single vortices cannot move away from these 
points even if one makes local gauge transformation, since if 
they did so, they would need an odd number of partners, as 
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Table II. The Z2 invariant numerically computed. The top two cases cor- 
respond to the uniform mass Wilson-Dirac model 8m = (b = 2 and t = 1 ) in 
the trivial C = phase, whereas the middle two cases are for C = ±1 phases. 
The bottom two cases are for the weak topological phase A and B . These six 
cases are indicated as black points in Fig. 2. 



m+ 


m_ 


C 


[ni»2»3«4] 


[[«in2«3»4]] 


-1 
9 


-1 
9 






[0110] 
[0110] 


[[0000]] 
[[0000]] 


0.5 

7 


1 

8 


1 
-1 


[1110] 
[0111] 


[[1000]] 
[[0001]] 


2.5 
3.5 


5.5 

7 






[1111] 
[0101] 


[[1001]] 
[[0011]] 



discussed above. Indeed, in the phase C - ±1 in Fig. 2, odd 
number of vortices are located on these four symmetry points 
{X/} in any gauge, as long as the phase convention (10) is used. 
See Table II. 

This implies that unpaired vortices on these points can be 
used to reveal the topological properties of the present sys- 
tem. Indeed, it informs us of the "parity" of the Chern num- 
ber. Moreover, different configurations of vortices on {X/} im- 
ply topologically different phases, since the location of these 
vortices are gauge invariant. Therefore, it is natural to expect 
that the C - phase can be further distinguished by the ob- 
struction on the four symmetry points. Notice that the present 
Hamiltonian 'H has inversion symmetry. Therefore, the ob- 
struction on the four points {X,} is associated with the parity 
of the wavefunctions. With respect to the parity operator P, 
eq. (10) implies that we choose 



iK~k) = PiKk). 



(12) 



At the point X\ — (0, 0), for example, this means that t/f(0) = 
Pi/f(Q). Let us assume that among the two occupied wavefunc- 
tions of t//(Q), n wavefunctions have the parity -1 and others 
have +1. Then, the former are the obstructions of the gauge- 
fixing condition (12). When n is even (« = or 2), there 
appear even number of the vortices on X\, which can move 
away from this point via suitable gauge transformation. How- 
ever, if n is odd (n = 1), one vortex is forced to locate on 
X\. The Z2 invariant in this meaning can be extracted from 
det[if/^ (0)Pif/(0)] = ±1, where the case -1 is the Z2 obstruc- 
tion of the condition (12). 

This simple observation is extended to other symmetry 
points Xj. Here, it should be noted that the wavefunction \p{k) 
is not periodic with respect to k x because of the twist operator 
U(k x ), as discussed below (10). In particular, ifr(—k x + n, k y ) = 
U(n)t//(-k x ,k y ) = U{7i)P\j/(k x ,ky), where ky = or n, and 
the phase convention at k x - n/2 is thus modified. Note that 
P - <x 3 ® 1 2 and U(n) = H2 ® 0-3 and hence U{n)P - <x 3 ®o"3. 
Thus, the condition (12), taking this periodicity into account, 
is explicitly written as 



t/riXi) = PuKXd, i =!,•••, 4, 



(13) 



where P t = crj, <g> 1 2 (i - L 4) and P, = cr 3 ® cr 3 (i = 2, 3). This 
leads to the Z 2 obstruction 



Si = deti/ f J (Xi)Pnff(Xi), 



(14) 



which can take + 1 only. Now let us define a set of four num- 
bers [ni«2«3H4L where Si - (-1)"'. These numbers are gauge 



invariant modulo 2, indicating that n, = 1 (0) imply obstruc- 
tion (no obstruction) on X,. 

In Table II we show several examples of this invariant. No- 
tice that even in trivial C — phases on the m + - m_ line, 
the obtained invariant is [0110], which seems nontrivial at 
first sight. It is, however, due to the effect of the twist op- 
erator U{n). Roughly speaking, U{n) — I2 ® 0-3 implies that 
among four states, half are periodic and the other half are anti- 
periodic in k x , and therefore, the half-filled ground states in- 
cludes one periodic and one anti-periodic states. Therefore, 
we should define the relative Z 2 invariant as "excitations" 
for this background obstruction such that [[ni«2«3«4]] = 
[nin2«3«4] - [0110] mod 2. 

This is the bulk Z2 invariant that characterizes the weak 
topological phases in the present superlattice system. From 
Table II, we can interpret the phase A as [ [ 1 00 1 ] ] = [ [ 1 000] ] + 
[[0001]] and thus the C - 1 + (-1) = nature of this states 
is established. It is also obvious that the phase B, assigned 
[[001 1]], is distinguished from this phase. 

Let us finally mention how the Z2 indices n,- are related to 
the edge spectrum in the two different ribbon geometries as 
shown in Fig. 3. Define S t - (-1)"'. This is the relative Sj 
in (14) with respect to the background. Then, in the spirit of 
Ref. 37 we consider a product of two 5,'s [S- k l and Si 2 ] on 
a line k^ - k M (k y = 0, n, while k x - 0,n/2), and define 
the quantity n\ = Sii^ki- Again, it\ - +1. If ttj, - — 1, 
the edge spectrum in the ribbon geometry directed to the yu- 
axis becomes gapless at k M - k^, otherwise the spectrum is 
gapped. Here, this statement can be verified empirically using 
the explicit values of «, listed in Table II, while the proof of 
the statement involves calculation of the Berry curvature inte- 
grated along the loop k M - k^ (Wilson loop). We leave detailed 
description of this proof to a forthcoming publication. 
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